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Abstract 



The combinatorial problem of counting the black hole quantum states within the 

u ; 

Isolated Horizon framework in Loop Quantum Gravity is analyzed. A qualitative 
understanding of the origin of the band structure shown by the degeneracy spec- 
trum, which is responsible for the black hole entropy quantization, is reached. Even 

oo ; 

OO . when motivated by simple considerations, this picture allows to obtain analytical 



' expressions for the most relevant quantities associated to this effect. 

PACS numbers: 04.70.Dy, 04.60.Pp 



"*-! ' I. INTRODUCTION 



In 1974 S.W. Hawking 



in 



established that black holes behave like black bodies in the 
thermodynamical sense. This remarkable work provides a clear evidence that the similar- 
ity between the laws of black hole mechanics 0| and the ordinary laws of thermodynamics 
is much more than a mere mathematical analogy. This physical analogy is summarized 
in the Generalized Second Law j^j], which endows black holes with physical entropy in a 
pure thermodynamical sense. Any quantum gravity theory proposal has to provide the 
microscopic degrees of freedom that account for that entropy. 

Loop Quantum Gravity (LQG) offers a detailed description of the black hole 
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horizon quantum states. Black holes within LQG are treated in an effective way in the 
Isolated Horizon framework introduced by Ashtekar et al In this framework the hori- 
zon quantum degrees of freedom are described by a U{1) Chern-Simons gauge theory and 
fluctuate independently from the ones of the bulk, giving rise to the black hole entropy. 
At present, two inequivalent proposals s[ for characterizing the black hole degrees 
of freedom have received most of the attention. It is interesting that, within both of 
them, the problem of computing the black hole entropy can be reduced to a well defined 
combinatorial problem. 

To exactly solve this combinatorial problem is, however, a rather non trivial task and, in 
order to obtain analytical solutions, some approximations have to be made. In particular 
the large area approximation permits to perform an analytic counting of the black hole 
microstates Q, 12]. Using this approximation the theory reproduces the semiclassical 



proportionality relation between entropy and area and gives an additional logarithmic 
term with a —1/2 coefficient, 



where 7 is the Barbero-Immirzi parameter jl3| (a free real parameter in the theory) and 
7 C a numerical constant obtained from the counting. Fixing 7 to be equal to 7 C ensures 
consistency with the Bekenstein-Hawking entropy law for large areas. An important fact 
is that both definitions for the horizon states to be considered agree with (DO), the only 
difference being the value of j c . 

Alternatively, the complexity of the combinatorial problem can be overcome by telling 
a computer to make an exact counting by explicitly enumerating all states Though 
the exponentially growing number of states limits the counting to modest black hole 
sizes (a few hundred Planck areas), the results in this regime agree with the analytical 
computations in the large area limit. Even more, this direct computation reveals a richer 
behavior shown by the spectrum when avoiding any approximation. The most degenerate 
quantum configurations accumulate around certain evenly spaced values of area, with a 
much lower degeneracy in the regions between those values, thus giving rise to an effective 
"quasidiscrete" equidistant area spectrum, despite the fact that the area spectrum in LQG 
is not equidistant. Furthermore, this phenomenon is independent on the particular choice 
for the characterization of the horizon degrees of freedom. This result provides a contact 
point between LQG and the Bekenstein's conjecture and has important consequences 
for the physical properties of black holes, such as the entropy, which displays an effective 
discretization [I6fl, or Hawking radiation, that could carry some quantum imprints coming 
from the horizon structure at the Planck scale 17l |. 
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This "band structure" arising in the black hole area spectrum of LQG calls for a more 
intuitive explanation, unraveling the origin of this phenomenon from the theory. This is 
the main goal of the present paper. A recent work in this direction has been done by 
Sahlmann in [181 ]. where he gives some quantitative information about the black hole area 



spectrum. In this paper we will follow a rather different approach. Despite the complexity 
of the combinatorial problem, which makes a meticulous analysis unfeasible, the states 
can be properly handled by attending to a few properties that allow us to obtain the most 
relevant qualitative and quantitative information about the area spectrum, shedding some 
light on its behavior. In particular, this approach will help us to understand qualitatively 
the origin of the "band structure" and will also allow us to compute analytically the value 
of area corresponding to each "peak" of degeneracy. 



We have organized the rest of the paper as follows. In sections (JTTJ) to (|IV|) we review 
the previous works, paying special attention to the aspects related with our arguments, 



and establishing the notation we are going to use, while sections (|Vj) and (IVII) contain 
the main body of the present work. Section (JTTJ) is devoted to set up the combinatorial 
problem. In section (|IIII) some previous analytical results are presented. Section (ITVl) 
contains a summary of the computational results that showed the behavior that we are 
going to analyze. We present our qualitative picture and our quantitative analytical 
computations in section ((Vj). The main results are analyzed in section (IVII) . We finally 
conclude with an outlook in section (IVIII) . 



II. COUNTING AND LABELING CHOICES 

In the Isolated Horizon (IH) framework in LQG black holes are treated in an effective 
way, since they are introduced from the outset as an inner boundary of the spacetime 
manifold before the quantization procedure is carried out (see 0, ^| for details). Isolated 
Horizon boundary conditions are then imposed, which translate into quantum boundary 
conditions after the quantization procedure. The horizon states are described by a U(l) 
quantum Chern-Simons gauge theory, while gravitational degrees of freedom of the bulk 
are represented by spin networks, a set of edges with spin-like quantum numbers (j, m) 
(j G Z/2, m = {—j, —j + 1, j}) that intersect to each other at vertices. When an edge 
of the spin network pierces the horizon creating a puncture, it endows it with a "quantum" 
of area given by 

a(j) = 8tt 7 4vW+1) > (2) 

where j is the corresponding label of the edge, and with a quantum of curvature given 
by the label m (since the Isolated Horizon boundary conditions relate this label with the 
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U(l) Chern-Simons states on the horizon surface). Then, the quantum states of a black 
hole with area A must satisfy that the sum of the contribution to the area from each 
puncture equals the total horizon area, 

p 

A = 87r 7 4X)Vii0i + l) > (3) 
i=i 

where p is the number of punctures on the horizon. Also a condition coming from the 
fact that the horizon is spherical must be imposed. This is 



X>* = ' (4) 



which is called "projection constraint". The problem of counting the black hole microstates 
that account for its entropy is now reduced to a mathematically well defined combinatorial 
problem which can be stated as: 

How many different configurations of labels distributed over a set of distinguishable 1 
punctures are there, for all possible finite numbers of punctures, such that the constraints 
(jSJ) and (p0) are satisfied? 

There exists a certain ambiguity at this point, since there are two proposals concerning 
which labels have to be considered to account for all microscopic configurations. The issue 
of which is the proper counting is, however, beyond the scope of this paper, as the behavior 
that we want to analyze is obtained within both of them. 

The first of the two proposals was done by Domagala and Lewandowski in Q] and was 



complemented by Meissner in There, it is claimed that the horizon states are given 
by punctures carrying only the m; labels (as these are the labels related to the horizon 
states through the IH boundary conditions). The constraint ([3]) is then reinterpreted in 
terms of ImA. 



The second proposal is due to Ghosh and Mitra [10] , and it considers that both labels, 
ji and ?7i j, characterize the horizon quantum states. In this case, both constraints (j3J) and 
Q can be imposed as written above. The structure, results and main differences between 
both models can be seen in Table fl] 

For the purpose of this paper, we need to deal with the labels related to area, so we 
will call this labels generically Sj, in such a way that the s, will correspond to \rrii\ in the 
first case and to ji in the second one. Furthermore, for the sake of simplicity, we will deal 
only with integer numbers, so that we will take Sj = 2\rrii\ or Sj = 2ji in each case. Then, 
in the DLM case, there will be two possible values of for each Sj, namely {— 



1 The fact that punctures are distinguishable has its origin in some subtleties related with the action of 
diffeomorphisms during the quantization procedure, and plays a key role in the combinatorial problem. 
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Table I: Comparison between the DLM and GM countings 





DLM 


GM 


Labels 


mi 




Area 




8vr 7 ^ E* ViiC/i + 1) 


Projection constraint 


Ei m i = 


Ei m i = 


Entropy 


5(A) = 7DMiA_l lnA 


5(A) = 2£M|-|lnA 


BI parameter 


7 = 7dlm = 0.23753 


7 = 7GM = 0.27407 



while in the GM one the possible values of m ; will be {— f-, — 4f + 1, 4f}, so there will 
be Sj + 1 values of m,; for each Sj. This will be the only difference that we will have to 
introduce in our arguments in order to account for both counting models. 



III. PREVIOUS ANALYTICAL RESULTS 

In this section we review briefly the previous analytical results on the counting of 
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black hole microstates ft la lid, I 111 . Il2| present in the literature. When addressing the 
combinatorial problem described in the previous section, a key point is to consider that 
punctures are distinguishable, as shown in fl. With this in tnind , one should consider 
all possible orderings of labels over punctures. But given a certain configuration of Sj 
labels, all possible reorderings give rise to states with exactly the same area. One can 
then characterize a configuration just by fixing the number n s of punctures that take 
each particular value of s and introducing all possible orderings as a certain degeneracy 
associated with this configuration. Thus, in the remainder of the paper, a given set of 
numbers {n s }^"^ (where s max is the maximum value of s) will be called configuration. 
A configuration will be permissible if it satisfies the constraint (j3J), which in terms of n s 
reads 

k 

^£lJ2 n sV^s + 2) = A • (5) 

s=l 

Then, in order to consider all quantum states contained in a given configuration, one has 
to take into account the degeneracy coming from two sources: 

• one due to all possible reorderings of the {s^ labels over punctures, 

• and the other coming from all possible combinations of the m ; labels associated to 
each configuration satisfying the constraint (TjJ . 
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The difference between the two possible countings is contained in this last term. For 
the only reason of being able to explicitly write down some expressions, we are going to 
consider for the moment the term corresponding to the GM counting. One can then write 
the degeneracy associated to a given configuration {n s } s s ™l x as: 2 

^i,.,0 = ^rn( s + 1 ) ns > ( 6 ) 

where sums and products run from s — 1 to s max . In the above expression the projection 
constraint is not being introduced, but this fact will not affect the results that we are 
going to obtain in the remainder of the section. This degeneracy was studied in 

HQ, 

where the question of which are the values of n s that give rise to the maximal value of 
degeneracy, for a fixed value of area, was addressed. 

The degeneracy d(nx, n Smax ) (or equivalently lnd(ni, n Smaa )) is maximized by 
varying n s subject to the constraint (JHJ), which is introduced via a Lagrange multiplier. 
This maximizing process can be easily worked out in the large area limit, where the vari- 
ables n s 3> 1 can be considered as continuous. The variational problem is then easily 
solved by using Stirling's approximation, which gives the result 3 

n* 



where for consistency, A must satisfy the "normalization condition" 

J2(s + l) e -V^) = lt (g) 

s=l 

Numerical solutions of this equation, in the large area limit (s max 3> 1) gives Aqm = 
0.861006 ( or Adlm = 0.746232 in the case of the other counting proposal). 

We will call this n s distribution the Maximal Degeneracy Distribution 4 ' (MDD), and it 
will play a pivotal role from now on. Besides, it was shown in HQ that the introduction 
of the projection constraint does not modify this distribution, so despite one starts without 
imposing it, the results can be considered as including this constraint. 

It is worth to note that, in the MDD, the proportions between the different n s are 
maintained for different values of area (the values of n s grow proportionally), and then 



2 The factor (s + is the one accounting for all the possible values of associated to each s,, so in 
order to make the analysis for the counting of [9l|_il will be enough to change this term by 2 n " 

3 Although the equivalent expression obtained in [llj was presented on the basis of different considera- 
tions, it can also be given the same interpretation of a degeneracy maximizing distribution. 

4 We will use the term distribution as opposed to configuration, in the sense that it gives the proportions 
between the different n s instead of the absolute value of each n s . We will use this terminology in the 
next sections. 
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the values of h s are independent of area. When plotting h s (Figured]), an interesting 
behavior is observed. Although the largest contribution comes from the smallest value of 
s (which contributes with approximately one half of the punctures) , the contribution of 
the next few values of s is also significant. Nevertheless, the MDD shows an exponential 
decrease as s grows, so for s larger than the smallest few values the contribution will 
become negligible. 
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Figure 1: The h s given by the Maximal Degeneracy Distribution (MDD) is plotted as a function 
of s. The relevant contribution of the lower values of s and the exponential decrease as s grows 
are observed. 



Once the MDD has been obtained, the tota 



value of area can be computed. The result is 



HQ 



number of quantum states for a given 



d= " e ^ , (9) 



where a ~ 0(1). It is seen that the number of quantum states grows exponentially with 
area; the extra factor A~ x l 2 appears due to the introduction of the projection constraint 
((4|). From this the entropy can be computed, obtaining 



s <^5H ln l +o( ^ 0) ' (io) 

This result verifies the semiclassical Bekenstein-Hawking entropy formula for large ar- 
eas provided that 7 = X/tt. Substituting the value of A for each counting the corresponding 
values for the Barbero-Immirzi parameter are obtained 

7gm = 0.274066858 , 7 DLM = 0.237532958 . 
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A. Large area limit 



In the previous computations the large area approximation was involved. However, 
one can wonder about the meaning of "large area" in this context. If one looks at the 
normalization condition ((HI), it is easy to see that the value of A obtained from it depends 
on the value of s max to which we are summing up. Then, as the value of s max depends on 
the area, we have a A that is a function of area. Nevertheless, if one studies the function 
A = X(A) (or equivalently A = A(s max ), as shown in Figure[2]), one sees that the value of A 
grows very quickly and saturates the asymptotic value for relatively small values of s max 
(for s max around 12 the value of A only differs from the asymptotic value in a 0.006%). 
But this value of s max corresponds to values of area around 45£p. So for areas larger than 
that, we can say that we are already in the large area limit, as far as the distribution (0) 
is concerned. 

1.0 - 
0.9 - 

• * * * * * * * " * * * * * * * 

0.8 - . 

A(S,„„) 0.7 : 

0.6 - 
0.5 - 

0.4 t ! i , i i i , i , i , i , i i i , i , L 

5 10 15 20 

Smax 

Figure 2: The value of A as a function of s max is plotted, and compared with the asymptotic 
value. 



IV. PREVIOUS COMPUTATIONAL RESULTS 

In the previous section, some approximations were employed in order to count the 
number of quantum microstates compatible with a macroscopic black hole. One can 
legitimately be worried about the fact that this approximations could be hiding part of 
the richness of the problem. Fortunately, in spite of its intrinsic complexity, an exact 
counting can be performed to see whether there is a richer structure in the spectrum or 
not. This can be done by means of an explicit enumeration computational algorithm. The 
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strategy is to generate systematically all possible combinations of labels (for any possible 
number of punctures), and to check one by one whether it satisfies the required conditions. 
Then, by explicitly enumerating all states, one can make an exact counting of the black 
hole quantum configurations (for a given value of area) in this framework. This was done 
in [141 ]; here we are going to review the main results obtained there and in subsequent 
work. Even when such an exact counting can be done, the price to pay for overcoming the 
complexity of the problem with an explicit enumeration is a severe restriction to the black 
hole sizes that can be analyzed due to the huge number of configurations to be counted. 
For that reason, the available computing power allowed to analyze black holes up to just a 
few hundred Planck area sizes. However, these computations were enough to confirm the 
results of the previous sections, namely the exponential growth of the number of states 
with area and, when imposing the projection constraint, the factor A~ 1 / 2 . The fact that 
this results are compatible with the analytical computations gives one some confidence in 
the interest of performing such a counting even though, due to computational limitations, 
one is restricted to work in a small horizon area regime, far below the large area limit in 
which the Isolated Horizon framework in LQG was originally formulated. 
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Figure 3: Plot of the degeneracy (number of different horizon states in each area interval of 
0.01£p). States accumulate around some equidistant values of area, exhibiting a band structure. 

But besides confirming the previous analytical results, the exact counting showed a 
much richer behavior in the black hole area spectrum [16||. It was found that the black 
hole quantum states are distributed according to a "band structure" in terms of the area. 
The most degenerate configurations cluster around evenly spaced values of area, giving 
rise to equidistant peaks of degeneracy, with some orders of magnitude less degeneracy 
in the regions between them (Figure 02). This fact gives rise to an effective equidistant 
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quantization of the black hole area in LQG, even when the area spectrum in the theory 
(j2J) is not equidistantly quantized. The most relevant quantitative information about this 
phenomenon is the fundamental area gap between peaks, which is given by 



AA 



7X4 , 



(11) 



where x was estimated to be 



V 



.80 



A remarkable fact is that this result was obtained for both choices of labels to be counted, 
and that all the difference resides just in the value of the Barbero-Immirzi parameter. 

The obvious interest is now in the physical consequences of this structure. The first 
clear consequence is in the entropy-area relation. This periodic band structure in the 
area spectrum gives rise to a very distinctive signal in the black hole entropy, namely a 
stair-like behavior of entropy function of shown in Figure [H Furthermore, 

the particular structure of the area spectrum can also have some implications regarding 
the black hole radiation spectrum, as pointed out in [17|. 
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Figure 4: Plot of the results for the entropy as a function of area (in Planck units) obtained with 
the computational counting. The stair-like behavior, with a step width corresponding to AA, is 
observed. 

On the other hand this regular pattern in the black hole area spectrum provides a 
nice contact point with the heuristic ideas of Bekenstein and Mukhanov [3| about black 
hole area equidistant quantization. Even though the basic area spectrum in LQG is 

5 For details on how to obtain the entropy shown in Fig|4]from the degeneracy of Fig[3]see 0, 0| 
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not equidistant, this phenomenon shows that in the case of black holes this equidistance 
in the spectrum appears in a rather subtle way, namely as a result of the non trivial 
degeneracy distribution. This point of contact becomes even more intriguing when one 
realizes that the value of x is close to 8 In 3 ~ 8.788898, as there is also a logarithmic 
constant arising from the heuristic considerations of Bekenstein and Mukhanov. It is 
evident that no reliable conclusion can be extracted from this numerical proximity but it 
is worth keeping it in mind to see if a more detailed work can confirm a deeper relation 
behind this coincidence. 

V. THE RICHNESS OF DISCRETENESS 

In this section we seek to understand where the equidistant structure in the black hole 
spectrum comes from. We are going to analyze what happens to the MDD when one 
takes into account the discrete nature of the problem. Then, we are going to classify all 
configurations in sets characterized by two parameters, in such a way that the accumula- 
tion of states around the peaks of degeneracy becomes explicit and easy to study. Using 
these parameters, and some information extracted from the MDD, we will compute the 
value of area corresponding to each peak of degeneracy and then the area gap between 
peaks. 

The first thing to consider is how is it possible to obtain information about the quasi- 
discrete structure of the spectrum using a distribution that was computed with approx- 
imations that seem to neglect all the information about this behavior. In this point the 
important thing to notice is that, in fact, the approximation that is hiding all the discrete 
information is to assume that one can find some configuration satisfying the MDD for any 
given value of area. When doing so, one is implicitly assuming that the n s numbers can 
take any possible value given by 0, that in general are not integer values. It is clear that 
a non integer value for n s makes no sense. Then in order to find the actual maximally 
degenerate configuration, one should take the closest integer to each value of n s given by 
the MDD. However, if one modifies the value of n s , one is modifying also its contribution 
to the area. Then, there are two possible cases, depending on the value of area we start 
with: 

• When one tries to find the closest integer configuration, the area changes of each n s 
compensate each other in such a way that at the end the integer configuration that 
we find takes almost the same value of area. Then we will be able to find some highly 
degenerate integer configurations with the same value of area we started from. This 
case would correspond with a peak of degeneracy. 
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• When changing to the closest integer configuration, the deviations of each n s add up 
giving rise to a global area change so that the resulting integer configuration lies, in 
fact, in a different region of the spectrum. If one tries to find some integer configu- 
ration with the same value of area than the continuous one, it will not suffice to just 
take the closest integer to each n s . One would be forced to modify considerably the 
n s distribution, in order to reach this value of area with integer n s values. But then, 
the obtained configuration will follow a distribution no longer close to the maximal 
degeneracy one, and would then have a much lower degeneracy. Therefore, one will 
not be able to find a highly degenerate integer configuration for this value of area. 
Such values of area are the ones corresponding to the regions of low degeneracy 
between peaks. 

Our task now is to find out which values of area correspond to the first case and which 
ones to the second. In order to do that, we will use a convenient classification of states. 

A. Classifying states 

The combinatorial problem we are trying to address is a very complicated one, given 
the large number of variables (degrees of freedom) that come into play. For this reason, it 
is very difficult to handle all the information in a straightforward way. In order to be able 
to understand the underlying structure, we are going to organize all these configurations 
according to two parameters that will allow us to have a reasonable number of variables 
while keeping enough information for our analysis and computations. The two parameters 
we are going to consider to classify configurations are: 

• The number p of punctures of the configuration, p = ^^=1* n s ? an d 

• the sum over all punctures of the s; labels S = Y^=i s i = J2l=T s n s- 

For each pair of values of these parameters, we will have a set of many possible con- 
figurations. But the interesting thing is that if one fixes a given pair (S,p), then the only 
freedom left to change the value of area associated to a configuration, is to distribute the 
S "units of s label" over the p punctures in different ways (or in other words, to change 
the n s distribution). But the changes in area given by changing the distribution of n s 
are very small compared with the change in area given by modifying the parameters S 
or p in one unit (the requirement that all n s must be integer obviously implies that S 
and p can only take integer values). Then, by considering all possible n s distributions, 
one is covering an almost continuous region of area in the spectrum, while modifying S 
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or p results in a discrete jump to another area region. Of course, if one modifies radically 
the n s distribution, from one extreme to the other, one can get changes in area larger 
than the one given by a change of one unit in S or p, so these different area regions could 
overlap at some points. 

On the other hand, although changing S produces a jump in areas and so does a 
change in p, one could in principle modify both parameters in such a way that the final 
area change is small. In fact, as we are going to see,there is a way of changing S and p so 
that the area does not change. As pointed out in [171], there is a very precise relation in 
the area spectrum of LQG that will help us to obtain this interesting relation between S 
and p. One can check that the contribution to area given by one puncture with s = 6 is 
exactly the same as the contribution given by four punctures with s = 1. The interesting 
fact about this relation is that it is the only existing one for the low values of s that are 
relevant to the highly degenerate configurations 6 (as pointed out in section IHIl the value 
of h s decreases exponentially with s in the MDD). 

Then, given a configuration, one can obtain another one with exactly the same value of 
area by removing a puncture with s = 6 and adding four punctures with s — 1 (decreasing 
the value of n 6 in 1 unit and increasing n\ in 4 units). But this change implies increasing 
the number of punctures p in three units and decreasing the sum of s over all punctures 
(S) in two. Therefore, different pairs of parameters (S,p) related by this transformation 
will be in the same area region. 

We can write down this relation in a more concrete way. Given a value So for S and a 
value 7 po = 1, 2, 3 for p, all pairs (S t ,pt) that satisfy the following relation: 



with t 6 Z such that S t > pt, are in the same region of area. Sq will be the maximum 
value of S and po the minimum value of p among all pairs (S t ,pt) satisfying this relation. 
Thus, if we consider the quantity K = 3S + 2p 0l we can associate to the same value of 
area all pairs of parameters satisfying 



Then, for each value of K we will obtain the configurations that appear in a certain 
region of area. In fact, it is important to notice that, if one takes into account the 



6 The next exact relation is found between one puncture with s = 16 and six punctures with s = 2, 
but the contribution of punctures with s = 16 to the highly degenerate configurations is completely 
negligible. 

7 Any value of p larger than 3 would be in correspondence with one of these tree values of po, i.e., a pair 
(S,p = 4) will correspond to (So = S + 2,po = 1), and so on. 



(S t ,pt) = (S - 2t,p + 3t) , 



(12) 



3S + 2p = K . 



(13) 
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projection constraint, then the value of S can only be even (for an odd value of S would 
imply that J2i m i can on ly ta ^e half-integer values, and then it could not be zero). With 
this in mind, K will only be allowed to take even values. 

B. Highly degenerate integer configurations 

Now, in order to account for the peaks of degeneracy, we need to consider the most 
degenerate integer configurations, which we will find with the help of the MDD. For a 
certain value of area, this distribution fixes a value for S and p, that we will call S m d(A) and 
p md (A). Furthermore, as the values of n s are constant with area, S m d(A) and p m d(A) will 
grow proportionally with area, giving rise to a constant quotient s m d = r 2 ^. Hence, not 

Pmd 

all the pairs (S,p) can contain maximally degenerate configurations; these configurations 
can only take values of S and p satisfying the quotient s md . However, it is important 
to notice that the values of S md (A) and p m d(A) fixed by the MDD are not, in general, 
integer numbers. Then, if starting from a configuration satisfying the MDD one changes 
to the closest integer values for each n s in order to find the actual maximally degenerate 
integer configuration, this would necessarily imply a change in S and in p to integer 
values. But as we have seen, to change S and p implies relatively large changes in area, 
unless such changes in S and p follow the "constant area" relation (fT3l) . Therefore, if those 
(S m d(A) , p m d(A)) satisfy this relation for an even value of K, it will be possible to find, for 
the same value of area, some integer pair (S',p'), with even S', close to (S m d(A) , p md (A)) , 
thus containing highly degenerate configurations. Otherwise, it will not be possible to find 
any highly degenerate configuration for that value of explained at the begining 

of section [V] 

In addition, among all configurations compatible with a given pair of values (S',p'), 
the most degenerate ones will be those having n s distributions close to the MDD, and 
therefore they will all appear together in a region of area much smaller than the total area 
covered by the set of all configurations with these values of (S',p'). Then, although the 
regions of area corresponding to different pairs (S, p) can overlap (as pointed out above), 
the regions containing highly degenerate integer configurations will not. Thus, we expect 
to find the highly degenerate configurations clustered around some area values, each one 
corresponding to a different value of K. 

We are going to compute in next section the values of area for which the MDD fixes a 
pair of values (S md (A) , p md (A)) satisfying the constant area relation for each value of K, 
which will correspond to the values of area of the peaks of degeneracy. 

We can understand the above discussion in a more graphical way looking at Figure 
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1 2 3 4 5 6 7 8 

V 

Figure 5: Plot of the sum S of spin labels vs. the number p of punctures. All the discrete 
configurations are placed in the marked points. The thin lines represent "constant area surfaces" 
(K-lines) while the thick line cointains the values of (S,p) that satisfy the quotient s m d (MD-line). 

El The positive slope line represents the pairs (S,p) that satisfy the maximal degeneracy 
quotient s m d (MD-line). Each of the negative slope lines represent the values (S,p) related 
by ( EDS ) for each even value of K (if -lines). The marked points represent the allowed 
integer pairs (S,p). Only for those points at where the MD-line intersects some of the 
-fT-lines, a close integer pair (S',p') can be reached following the corresponding K-\me 
(keeping area constant). In other words, the actual highest degeneracy configurations, 
that can only be found in the integer points close to the MD-line, correspond to the 
values of area associated to the point at which the corresponding K-\ine intersects this 
MD-line. The MDD provides the necesary information to compute the value of area 
associated to each point of the MD-line, as well as the slope of this MD-line. We can 
compute, therefore, the area of each intersection point and, thus, of the corresponding 
peak of degeneracy. Finally, computing the difference in area between two consecutive 
intersections (two consecutive even values of K) we will get the area gap between two 
peaks of degeneracy. 

As a final remark in this section, we can analyze the effect of the projection constraint. 
In our model, this constraint is introduced by considering only even values of K (i.e. even 
values of S). Then, if the projection constraint was not introduced, there would be an 
additional line between each two consecutive if-lines in Figure [5l This would correspond 
to having an additional peak of degeneracy between each two. But looking at Figure 
[3] one can see that, given the proportions of the spacing between peaks and the width 



15 



of those peaks, placing an additional one between each two would almost result in no 
low degeneracy regions between them, hiding then the "quasi-discrete" structure of the 
spectrum. Then, as pointed out in [I8fl, the regular pattern we are studying is a general 
feature that affects to all states and not only those satisfying the projection constraint. 
But it is precisely the introduction of this constraint what makes the "discrete" structure 
to arise in a clear and relevant way. 



C. Computation of AA 

Let us then proceed to the explicit computation of these quantities. The steps we are 
going to follow are: 

• In the first place, using the MDD we will compute the quotient for maximally 

Pmd 

degenerate states (s m d). 

• From ( |T3l) we will obtain an explicit relation S = S(p, K). 

• We will use this explicit relation and the value of the quotient s m d to compute the 
number of punctures of the maximally degenerate state p m d(K) for a given value of 
K (the value of p at which the MD-line intersects a K-\me in Figure [5]) . 

• Then, again using the MMD, we will compute the mean contribution to area A m d 
of a puncture in a configuration satisfying this distribution. 

• Thus, the value of area associated to an intersection with a line characterized by K 
will be A md (K) = p md {K)A md . 

• Finally, computing the difference between A md (K) for two consecutive even values 
of K (A md (K + 2) - A md {K)), we will obtain the value of AA. 

In order to compute s m d it is worth noticing that the quantity ^ can be seen as the 
mean value of s of each puncture in a configuration. Then, to compute this value in the 
case of the MMD we can write 

Smd = ^2 sn s . (14) 

s 

Thus we can compute the value of s m d and we know that 

= Smd ■ (15) 

Pmd 
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Now, from the relation between S and p in a given band (3S + 2p = K), we can extract 
the following equation 

S( P ,K) = ^-^p. (16) 



Plugging this into ( jl5h . we get 



Pmd(K) Pmd(K) 



leading to 



We have then the number of punctures that correspond to a maximal degeneracy config- 
uration for a given value of K (the intersection for a given i^-line) . 

Now, to compute the mean contribution to the area from a puncture in a configuration 
satisfying the MDD, we proceed in the same way as we did to compute s md . We then 
write 

And = = 47F ^P £ ™ s V S ( S + 2 ) ■ (18) 

s s 

With this expression we can write the value of area associated to each of the intersections 
for each value of K, 

A md {K) = Pmd(K)Amd = m " ■ (19) 

We have then arrived at the main goal of the paper, i.e. obtaining the value of area 
associated to the corresponding peak of degeneracy for each value of K. With this expres- 
sion we can compute numerically the value of area of each peak of degeneracy. We can 
also easily see that A md has a linear dependence on K, so the peaks are evenly spaced. 
We can hence compute this spacing just by taking the difference between two consecutive 
values of K, 

2A 

AA = A md (K + 2) - A md {K) = A md (p md (K + 2) - p md {K)) = md . (20) 

<jS md + ^ 

Then, finally, writing explicitly all the terms in the above result, we can express the value 
of the area gap between peaks as 

a a SttWI Es V^+2)(* + l)e- A ^ 

AA GM = Xgm7gm = t== • (21) 

3(E s s(s + l)e- AGM V^+ 2 )) + 2 

At this point, we can recall that the only difference in all this discussion between the 

label choice we are using and the other comes from the degeneracy associated to the 
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combinations of m, compatible with each configuration. Introducing this change, the 
result for the AAdlm in the case of the counting of {s| is 



AA v ^ SttWI E 5 2 ^g+ge^V^ 

^^iDLM — XDLM7DLM — , = , (/^J 

with the corresponding Adlm- 



VI. ANALYSIS OF THE RESULTS 

In this section we present the numerical values obtained for % aR d analyze them. 
The resulting values of the expressions we have found can be easily computed using 
Mathematica™ and we get 

Xgm = 8.789242 , X dlm = 8.784286 . 

The fact that the difference between these two values is in the fourth digit gives us a hint 
on the level of accuracy that is being reached. Besides, it was pointed out in {3| that the 
value of x is numerically close to 8 In 3 = 8.788898. One can see that the above results 
coincide, also up to the fourth digit, with this value, and furthermore, that the value of 
8 In 3 is contained between the two above values of x- One can compute the deviations 
between those three values: 



|xgm - 81n3| 

81n3 
Ixdlm - 81n3| 

81n3 
Ixgm — Xdlm| 



0.000039 = 0.004% 
= 0.00052 = 0.05% , 
= 0.00056 = 0.06% . 



Xgm 

Then, with a precision of 0.06%, the values of Xgm, Xdl and 8 In 3 are the same. Of 
course, this is still not a rigorous proof that x is equal to 8 In 3, but it is relevant to see 
how, when one improves the accuracy of the calculations, the numerical coincidence keeps 
being satisfied. 

Let us end this section with two remarks. 

• One can check whether the results obtained with the model presented here are in 
good agreement with the computational data obtained from the algorithm of 
It can be seen in figure ([6|) how the values for the area of the peaks that we obtained 
here fit the peaks observed in the spectrum obtained from the computer. We see 
how the analytical values match in a nice way with the computational data. On 
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Figure 6: Comparison between the analytical values for the area of the peaks (dashed vertical 
lines) and the actual peaks obtained from the computational data. 
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Figure 7: The mean values of n s obtained from the computational results (points) are compared 
with the exponential decrease expected from the MDD. 

the other hand, in figure (J7]) the mean values of the h s obtained with the computer 
for the five most degenerate configurations of three consecutive peaks (with areas 
between 170 and 177£p) are compared with those given by the MDD. One can 
check that, even for this extreme low value of area, the agreement is quite good, 
as expected from the analysis in section IIII Al so one can feel confident with the 
use of the MDD in the computations. Finally, using these computational data we 
have observed that, in fact, all configurations giving relevant contributions to the 
degeneracy at any given peak are characterized by pairs of values (S, p) that satisfy 
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the relation (TT3]) for the corresponding value of K, in complete agreement with the 
analysis presented in section [V] Thus, the computational data support the fact that 
the model presented here works reasonably well. 

• At this point we can analyze the results previously obtained in There, the 

problem is addressed using a rather different approach, namely, reformulating it 
in terms of the so called random walks. It is very interesting to see how, within 
this alternative approach, the accumulation of states around certain values of area 
also becomes manifest. By treating the area spectrum of LQG as an effectively 
quasi-equidistant spectrum, a way to compute the area gap between peaks is then 
proposed. The value for A A was obtained as 2/3 of the spacing in this effectively 
quasi-equidistant area spectrum. This 2/3 factor was introduced ad hoc in order to 
fit the computational data. A noteworthy fact is that this independent derivation 
gave rise to the same expresion (1221 . Nevertheless, from the point of view of the 
authors it is not easy to reconcile the introduction of this 2/3 coefficient with the 
qualitative picture of a quasi-equidistant spectrum being the origin of the observed 
regular pattern. 

With the picture presented here, it is now rather easy indeed to understand where 
this coefficient comes from. The fundamental area gap in the quasi-equidistant 
spectrum of |3| corresponds to the mean area change given by increasing S in 
one unit in our formalism. But as we have seen, S only takes even values. Then 
the minimum increment in S must be two. Furthermore, if one comes back to the 
relation between So, po and K (K = 2So + 3po) then one can check that for each 
even value of Sq there are three corresponding values of K (the ones corresponding 
to po = 1,2,3) and then to three peaks of degeneracy. Hence, the mean area 
change given by increasing S in two units corresponds to three times the area gap 
aetween peaks. Thus, the fundamental gap of the quasi-equidistant spectrum in 
is nothing but three halves of the area gap between peaks A A. 



VII. CONCLUSION AND OUTLOOK 



Let us summarize the results of the paper. We have analyzed the combinatorial problem 
and we have qualitatively understood the reason why the highest degenerate configurations 
can only appear for some values of area and not for all of them. When the discrete nature 
of the problem is taken into account, there are regions of area for which the "discrete 
configurations" are not allowed to satisfy a distribution close to the one that gives the 
maximal degeneracy in the continuous case, thus giving rise to the observed pattern in 
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the black hole area spectrum. We have also verified that the analysis is valid for both 
choices of labels, as the arguments presented here apply equally to both cases, so it seems 
now rather natural that the analyzed behavior of the spectrum appear with both counting 
procedures. Finally, our analytical computations allowed us to obtain the values of area 
for which the peaks of degeneracy should appear and showed that these values are evenly 
spaced. In addition, the results match in a nice way with the computational data obtained 
in [l^j], thus indicating the validity of the model. From this, we have also been able to 
compute the analytical value of the corresponding parameter \ for both label choices and 
we have found that the results coincide up to a precision of 0.06%. Furthermore, we have 
checked out that, up to this improved precision, the surprising numerical coincidence with 
8 In 3 keeps holding. 

There are still some important open questions. On the one hand, one may ask which 
are the sources of this 0.06% deviation. Moreover, it would be very interesting to obtain 
an analytical proof for the conjectured value of x — 8 In 3. On the other hand, although 
the area gap between peaks obtained with our model has no dependence on the area, one 
may be interested in knowing what would happen to the width of the bands for large 
areas. Whether this width increases with area, thus hiding the quasi-discrete behavior, 
or not, is also an interesting issue to be investigated. A comprehensive analysis of the full 
combinatorial problem could shed light on some of these questions. But undoubtedly, the 
most important and interesting open question is to find a consistent physical interpretation 
to this intriguing behavior of quantum black holes. 
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